Coalescence of the Fermi-surface-related diffuse intensity peaks in disordered alloys 
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The possibility of disappearance of the diffuse-intensity peak splitting induced by the Fermi 
surface (i.e., of coalescence of the intensity maxima) with decreasing temperature is predicted. The 
underlying mechanism is the compensation of the reciprocal-space curvatures of the self-energy and 
the interaction. The theory also describes similar results obtained earlier for two low-dimensional 
models with competing interactions. The coalescence is compared with the recently observed "ther- 



mal" splitting in Pt-V which can be explained in the same way. 
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In the recent paper (Tsatskis 1998a) we proposed an explanation of the temperature dependence of the diffuse- 
intensity peak splitting found for the disordered CU3AU alloy (Reichert, Moss and Liang 1996; see also Reichert, 
Tsatskis and Moss 1997). The splitting is associated with the Fermi surface (FS) effects and is observed also for 
other FCC alloy systems: Cu-Al (Scattergood, Moss and Bever 1970, Schonfeld et al. 1996), Cu-Pd (Ohshima 
and Watanabe 1973, Ohshima, Watanabe and Harada 1976, Saha, Koga and Ohshima 1992), Cu-Pt (Ohshima and 
1 i Watanabe 1973, Saha and Ohshima 1993), etc. It is the consequence of the indirect interaction between alloy atoms 
via conduction electrons in those cases where the FS has flat or nested (i.e., identically curved) areas (e.g., Krivoglaz 
1 1996). The resulting effective long-range pairwise interatomic interaction is characterized by minima at some positions 
between the X = (110) and W = (lj0) special (Lifshitz) points (SPs); this feature is then reflected in the short-range 
order (SRO) diffuse intensity (Fig. |l|). According to the proposed theory, the temperature-dependent splitting is the 
non-mean-field effect and the result of the wavevector dependence of the so-called self-energy. The self-energy £ of 
1 . the pair correlation function (for detailed discussion and relation to thermodynamics, see Tsatskis 1998b) enters the 
expression for the SRO intensity, 

o. 

/(k) = c(l - c) [-S(k) + 2/3V(k)} ' (1) 



where k is the wavevector, /(k) is the intensity in Laue units, c is the concentration, (3 = 1/T, T is the temperature in 



energy units and V(k) is the Fourier transform of the pair ordering potential V i3 = (V£ A + V? B )/2 - V,j B . Potential 
V^j 3 corresponds to the interaction between an atom of type a at site i and an atom of type /3 at site j. 

The k-dependence of the self-energy leads to the temperature-dependent shift of the /(k) peak with respect to the 
corresponding minimum of the interaction V(k) and, as a result, to the temperature dependence of the splitting. The 
positions of the I(k) extrema are determined by the condition dkl = which gives 



2d k V = Td k T>. (2) 



Here k is the wavenumber along a line (e.g., (hlO) or (IkO) in Fig. [j]) containing off-SP peaks; in what follows only the 
I(k) profile along this line is analysed. Eq. (||) shows that the k-dependent self-energy shifts the intensity peak position 
away from the V(k) minimum. The right-hand side of Eq. (||) contains T as a factor, and S(k) itself is a function of 
T, while the left-hand side is T- independent. The shift and the splitting (which are linearly related) depend therefore 
on temperature. In the two widely used theories of SRO, the mean-field Krivoglaz-Clapp-Moss (KCM) approximation 
(Krivoglaz 1969, Clapp and Moss 1966, 1968) and the spherical model (SM) (e.g., Brout 1965), the self-energy is 
k- independent, the right-hand side of Eq. (Q) vanishes, the 1(h) peaks coincide with the minima of VQt), and the 
. splitting does not change with temperature. 

This treatment, however, implicitly assumes that in Eq. (|l|) the k-dependence of the interaction term 2/3 V(k) in 
the area of the splitting is dominant. In this case the profile of the intensity closely follows that of the interaction, 
and there exists one-to-one correspondence between the V(k) minima and the I(k) peaks. The k-dependence of the 
self-energy in this part of the reciprocal space is relatively weak, though qualitatively important for the description of 
the temperature-dependent splitting. Such assumption is certainly correct at sufficiently high temperatures, where the 
KCM approximation (in which the self-energy is k-independent) works reasonably well. Meanwhile, as temperature 
decreases, £(k) grows faster than 2/3V(k), since the first correction to the KCM self-energy is of order ((3V) 2 (Tsatskis 
1998a, b). We can then encounter a situation where the variations of S(k) and 2/3F(k) with the wavevector are 
comparable. With temperature decreasing further, the k-dependence of the self-energy may even become considerably 
stronger; the I(k) peak positions would then be determined by the S(k) maxima. 
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FIG. 1. The reciprocal-space geometry of the splitting. The (hkO) plane is shown. Filled circles mark the positions of the 
SRO diffuse-intensity peaks. Small open circles correspond to the SPs X and W. 

The self-energy and the interaction arc, in general, qualitatively different functions of the wavevector. In particular, 
there is no special reason to expect £(k) to have any extrema away from the SPs. Therefore, at lower temperatures 
I(k) might no longer exhibit features characteristic for V^(k). We assume that (i) E(k) has extrema only at the SPs, 
and (ii) its variation with k in the area of the splitting becomes more and more important in comparison with that 
of 2/3I/(k) as temperature decreases. Then the qualitative picture of the temperature behaviour of the splitting is 
as follows: At high temperatures S(k) is almost k-independent, and the /(k) peak positions deviate little from the 
minima of V^(k). As temperature decreases, the k-dependence of X becomes more pronounced; the peaks move farther 
away from the V(\s) minima and towards that SP krj (either X or W) at which S(k) has a maximum. Eventually, as 
temperature reaches certain value To, the intensity peaks coalesce at this SP and the splitting disappears (Fig. ||). This 
effect has never been observed experimentally, though Reichert et al. (1996) assumed that (i) the coalescence took 
place at the first-order transition temperature T t (denoted there by T ) and (ii) the splitting grew with increasing 
temperature as (T — T t ) s ; the bifurcation exponent s = 0.38 ± 0.15 was determined by fitting to experimental data. 
The coalescence temperature T can be found from the condition of vanishing second derivative of /(fc) at the SP fco, 
since its sign controls the presence or absence of the splitting. At the SPs all the first derivatives vanish, and from 
Eq. @) it follows that 

(dll) ko = c(l - c)/ 2 (fc ) (d 2 £ - WdtV) ko . (3) 

The splitting therefore disappears when curvatures of the self-energy and of the interaction term at the SP fco com- 
pensate each other: 

2(dlV) ko =T(dlE) ko . (4) 

To analyse qualitatively the behaviour of the splitting close to the coalescence point, it is convenient to use the 
analogy with the Landau theory of second-order phase transitions (e.g. , Landau and Lifshitz 1980) . In this temperature 
range the splitting above To is small, and expansions of V(k) and S(fc) in powers of the deviation 5k = k — fco of the 
wavenumber from the SP fco can be used. To describe the split minimum of ^(fc), only the second- and fourth-order 
terms are needed; since a SP serves as the origin, the expansions do not contain odd powers of 5k. It is therefore 
assumed that in the area of the splitting V{k) and S(fc) have the following approximate form, 

f(k) = f(0) + a f (5k) 2 /2 + b f (5k) 4 /4, (5) 

where / = V, S, ay < 0, by > 0, as < (S(fc) has a maximum at fc = fco), and the sign of &s is arbitrary. Substituting 
Eqs. (||) for V and S into Eq. (|l|), we get the same expression (||) for the inverse intensity G — I -1 , where at T = To 
the second-order coefficient — c(l — c)(— + 2f3ay) vanishes (see Eq. (0)), while the fourth-order one be remains 
positive. In the vicinity of Tq we can then write ac = og(T — To), clq < and regard be as temperature- independent. 
The inverse intensity I~ 1 (k) thus behaves in almost the same way as the Landau free energy, and 5k is the analogue 
of the order parameter. The only difference here is that the role of temperature is reversed: I~ 1 (k) has a double 
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FIG. 2. Schematic temperature dependence of the SRO intensity profile (a) and of the peak positions (b). 



minimum above the coalescence temperature To and a single one below it. From this analogy it follows that the 
splitting increases with temperature as (T — To) 1 / 2 at small positive values of T — Tq. Unlike its Landau-theory 
counterpart, the obtained bifurcation exponent s = 0.5 is exact, since the intensity is a regular function of the 
wavevector expandable into the Taylor series. At higher temperatures the behaviour of the splitting changes, and its 
value starts to approach that of the splitting in V(k). In this regime the difference between the two decreases with 
temperature as 1/T, unless the alloy is equiatomic, in which case the decrease is faster (Tsatskis 1998a). 

Returning to the quantitative description of the temperature dependence of the splitting in the whole range of 
temperatures, we note that it can be considerably simplified by taking account of the particular behaviour of V(k) 
and S(fc) (Tsatskis 1998a). The Fourier transform F(k) of an arbitrary real-space matrix F defined on the FCC 
lattice varies with k as follows: 1st coordination shell does not contribute to its fc-dependence; 2nd and higher shells 
lead to the term proportional to cos2-7rfc; starting from 8th shell, cos47rfc term appears; 21st shell produces cos67rfc 
term, etc. Here k is defined as the deviation from the SP X and measured in the reciprocal-lattice units (r.l.u.). As 
a result, under the assumption that matrix elements of F beyond 20th coordination shell are of no importance, 



F(k) = A F + 2B F cos2nk + 2C F cos Airk 



(6) 



where the coefficients are linear combinations of the matrix elements Fij ; for explicit expressions for the relevant 
quantities Bp and Cp, see Tsatskis (1998a). To find the I(k) peak positions, we insert Eqs. (||) for V and £ into 
Eq. (||) for the I(k) extrema. The presence of the off-SP minima in V(k) implies Cy > 0, \By\ < 4Cy, while the 
assumption of the absence of such extrema in E(fc) requires \By\ > 4|Cs|. The positions kj of the intensity peaks 
away from the SPs above Tq are given by 



cos 27rfcr 



(2B V -TB S ) /4(2CV - TCfe) , 



(7) 



whereas the V(k) minima ky can be obtained from Eq. ([?]) by putting B^ = Cs 
into Eq. (Q), we obtain a simpler equation for the coalescence temperature To, 

2 (±B V + AC V ) = T (±B S + 4Cs) Tn , 



0. Substituting the same Eqs. 



(8) 



where upper and lower signs correspond to the coalescence at the SPs X and W, respectively. Contrary to Eq. (0), 
the explicit expression for To cannot be obtained, since By: and Cs are temperature-dependent. 

To calculate By; and Cs, a particular approximation for SRO must be used. Here we turn to the simplest theory 
leading to the wavevector dependence of the self-energy, the high-temperature expansion (HTE) (Tsatskis 1998a, b). 
The HTE is the expansion in powers of f3V; the second-order HTE approximation for the self-energy gives 



S(k) = S d + 2xP 2 W(k) , 



(9) 



) 2 ,w. l3 = v 2 , 

not contribute to the k-dependence of S. Then 



where x = (1 — 2c) 2 , Wij = Vij, and is the diagonal part of the self-energy in the site representation which does 



B E = 2x(3 2 B 



C E = 2x(3 2 C w 



(10) 
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FIG. 3. Positions of the intensity peaks versus dimensionless temperature t = 4C'vT/\Bw\ for the A3B alloy with the 
interaction V(k) such that Cw = and kv = 0.25 (i.e., the minima of V(k) are located half-way between SPs X and W). For 
this set of parameters = ig = 0.25, and the character of the temperature dependence of the splitting is defined by the sign 
of Bw- (a) Bw > 0, case (i) in the text, (b) Bw < 0, case (ii). 



where \B\y\ > 4|CV|. Substitution of Eqs. (10) into Eq. (0) leads to the result 

cos27rfc/ = -{B v - xf3B w ) /A(C V - x(3C w ) ■ (11) 

Calculation of the temperature derivative shows that in this approximation the behaviour of the splitting is determined 
by the sign of the T-independent quantity By/Cy — ByCw- 

df 3 (cos2Trk I )^x{B w C v ~B v C w )/4(C v -xPC w ) 2 . (12) 

(i) BwC'v > ByCw- The derivative (jlj|) is positive, and cos27rfc/ increases with decreasing temperature; the /(k) 
peaks move away from the V(k) minima and towards the SP X (Fig. ||(a)). Finally, the peaks coalesce at X when 
cos2nki reaches the value +1 at the bifurcation point 

T<f = x (B w + 4C W ) I (By + 4CV) . (13) 

(ii) B\yCy < ByCw- The intensity peaks shift towards the SP W (Fig. ||(b)). The splitting disappears when 
cos27rfc/ = — 1, i.e., at the temperature 

T W = x (B w — AC W ) I (B v - 4C V ) • (14) 

(Hi) BwCy = ByCw, or the alloy is equiatomic (c = 0.5): The derivative ( |l2|) vanishes and the second-order 
HTE approximation predicts temperature-independent splitting. In this situation the temperature behaviour of the 
splitting is defined by the third- and higher-order terms in the HTE for the self-energy. 

The predicted coalescence of the FS-related intensity peaks can, in principle, be described by using any theory of 
SRO which leads to the k-dependent self-energy. As was already noted, the frequently used KCM and SM approxi- 
mations fail to reproduce this feature of the exact self-energy. Apart from the HTE, the available theories include the 
cluster variation method (CVM) (Kikuchi 1950, 1951), the 7-expansion method (GEM) (Tokar 1985, Tokar, Masan- 
skii and Grishchenko 1990, Masanskii, Tokar, and Grishchcnko 1991), and the a-expansion (AE) (Tsatskis 1998a) 
closely related to the GEM. The corresponding approximations for the self-energy were compared by Tsatskis (1998b). 
However, the CVM is a direct-space method and as such is hardly applicable to the case of long-range interactions. 
To describe the off-SP minima of V(k), at least eight first potentials are needed, and accurate sets of inverse Monte 
Carlo (MC) interactions determined from experimental data are yet significantly larger (Schonfeld et al. 1996). The 
HTE approach is limited to the same extent as the KCM approximation; it leads to quantitatively correct results 
only at sufficiently high temperatures, as follows from its very name. The GEM utilizes the assumption about rapid 
decrease of interatomic correlations with distance, which is invalid for long-range interactions. This leaves the AE, a 
modification of the GEM specifically designed to be used in such cases (Tsatskis 1998a). Similarly to the GEM, the 
AE is expected to be sufficiently accurate at realistic temperatures. 
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We now compare the effect of the coalescence with another anomaly in the SRO scattering, the "thermal" splitting 
of the intensity peak, observed for the Pt-V alloy system and reproduced in the MC simulations (Le Bolloc'h, Cren, 
Caudron and Finel 1997). For this system, the (100) peak splits along the (hOO) line with increasing Pt content. 
A theory of this anomaly was proposed by Tsatskis (1998c), who also showed that the splitting can occur when 
temperature decreases at fixed composition. This prediction was confirmed by the MC calculations (Le Bolloc'h 
1997), though the effect has not yet been observed experimentally. According to the proposed theory, the thermal 
splitting can be explained by the same mechanism of the compensation at the SP of the S(k) and 2/3V(k) curvatures. 
However, in this case the interaction V^k) has a simple minimum at the corresponding SP and, in agreement with 
the mean-field arguments, produces the I(k) peak with no fine structure at higher temperatures. This peak splits as 
temperature decreases and the positive curvature of E(k) exceeds that of 2f3V(k). In the event of the coalescence, on 
the other hand, the fine structure of V(k) is erased at lower temperatures due to the growing negative £(k) curvature. 
Both phenomena can be conveniently described within the same Landau-type approach outlined above. The only 
formal difference is the sign of da- «g > for the thermal splitting and da < for the coalescence. 

The present theory shows that the coalescence temperature To does not, in general, coincide with the transition 
temperature T t and is more or less arbitrary. Nevertheless, it can accidentally be close to T t (unless the ordering 
transition occurs first). This is a possible reason for the rather good agreement between the theoretical value s = 0.5 
of the bifurcation exponent and the experimental result s = 0.38 ± 0.15 obtained for the CU3AU alloy (Reichert et al. 
1996). The remaining discrepancy can be attributed to the difference between T and T t , and also to the data fitting 
over the wide temperature interval (about 140 K). 

Finally, we note that the proposed mechanism of the curvature compensation explains also the results obtained 
for the exactly solvable ID Ising model with nearest- and next-nearest-neighbour interactions (Kulik, Gratias and de 
Fontaine 1989) and for the 2D ANNNI model in the framework of the CVM (Finel and de Fontaine 1986). Observed 
for these two models is, in fact, the same effect of the coalescence as described here: the I(k) peaks produced by the 
off-SP minima of V(k) move towards a SP as temperature decreases and transform at this SP into a single peak at some 
temperature before the onset of long-range order (for the ID model, at non-zero temperature). The exact solution 
in the former and the CVM in the latter case lead to the k-dependent self-energy whose SP curvature compensates 
that of the interaction term at the coalescence temperature. In all three cases the real-space interactions produce 
minima off SPs in the reciprocal space; however, Finel and de Fontaine (1986) and Kulik et al. (1989) considered the 
short-range competing interactions which are very different from the long-range interactions associated with the FS 
effects. 

In summary, we have pointed out a possibility for the FS-related splitting of the SRO intensity peak to disappear at 
some temperature above the order-disorder transition. The off-SP intensity peaks reflecting corresponding minima of 
the effective pair interaction in the k-space can move towards one of the SPs (X or W) as temperature decreases and 
eventually coalesce at this SP, transforming themselves into a single peak. The driving force behind this effect is the 
k-dependence of the self-energy which becomes stronger with decreasing temperature. The coalescence occurs when 
the self-energy and interaction curvatures at the SP compensate each other. The shift of the intensity peaks with 
decreasing temperature towards X (an increase of the splitting with temperature) was observed for the CU3AU alloy 
(Reichert et al. 1996), while neither the shift towards W (a decrease of the splitting with temperature) predicted by 
Tsatskis (1998a) nor the coalescence at any of these SPs were found. The coalescence can be described qualitatively 
by the Landau-type theory and quantitatively by any theory of SRO which is sufficiently accurate and leads to the 
wavevector dependence of the self-energy. This excludes two popular approximations, the KCM and the SM, in which 
the self-energy is site-diagonal. The effect of the coalescence has been compared with the "thermal" splitting of the 
SRO intensity peak which was observed experimentally and in the MC simulations and can be described within the 
same theoretical framework. 
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